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C^ ■ We revisit the Anderson localization problem on Bethe lattices, putting in contact various 

aspects which have been previously only discussed separately. For the case of connectivity 
3 we compute by the cavity method the density of states and the evolution of the mobility 
|/-\ , edge with disorder. Furthermore, we show that below a certain critical value of the disorder 

*^ ■ the smallest eigenvalue remains delocalized and separated by all the others (localized) ones 

by a gap. We also study the evolution of the mobility edge at the center of the band with 
the connectivity, and discuss the large connectivity limit. 
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§1. Introduction 



The effect of quenched disorder on non-interacting electrons can be dramatic. 
C^ ' The wave-function completely changes form going from an extended plane wave at 

zero disorder to a localized wave at large disorder. Actually, even an infinitesimal dis- 
order is enough to induce localization of all eigenstates in one and two dimensions.'^ 
pj ' Localized electrons do not lead to electronic transport except if the system is cou- 

O . pled to a thermal bath, as a consequence conduction properties changes completely 

in presence of localization. 

This phenomenon, called Anderson localization, was discovered in 1958 by P.W. 
Anderson. 2' Since then a huge amount of work has been devoted to it, see for 
CN I instance the brief review "Fifty years of Anderson Localization"'^ and the mono- 

^' ■ graphs .'21' '12' '13 'IS|I Exactly solvable models naturally played an important role in its 

f— ^ . understanding. One dimensional systems have been studied thoroughly.^' However, 

in these systems all states are localized and hence some aspects of Anderson localiza- 
("^ • tion, such as the transition between localized and extended states, and the critical 

O . properties at the mobility edge could not be studied. All these properties can instead 

be analyzed on Bethe lattices (see the next section for a precise definition), which 
therefore provide a very useful benchmark since they are still simple enough to be 
^ . studied without making any approximation. 

^ I The first analysis of Anderson localization on Bethe lattices was performed by 

Abou Chacra, Anderson and Thouless^''^ and many other studies followed both in 
the physics (see for instance the work3^'121''lIIJ.lII|i ^nd references therein) and math- 
ematic3l2J>llIli.ll3>ll5J communities. Our main goal in this work is to combine various 
point of views on the problem, and to underline a particular feature of Bethe lattices. 
In fact, besides the properties of the density of states, the evolution of the mobility 
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edge with disorder and with the connectivity k + 1, we shall discuss the presence of 
isolated eigenvalues. Without disorder, the smallest eigenvalue for the tight-binding 
Hamiltonian (or the discrete Laplacian) on a Bethe lattice is isolated from all the 
others by a gap.l^ This eigenvalue plays a very important role in determining the 
physical properties of some statistical mechanics models displaying a condensation 
transition at low temperature, such as the spherical model^ or non-interacting 
bosons on a Bethe lattice. The evolution of the isolated eigenvalue in presence of 
disorder is interesting in particular in connection with studies of disordered bosonic 
models since it may lead to somewhat different behaviors than the ones obtained 
for finite dimensional systems, where no isolated eigenvalue is found (even though 
the relevance of the spectrum of the one-particle kinetic energy is less obvious in the 
presence of interactions). 

§2. Model and iterative equations 

The Anderson model corresponds to the Hamiltonian: 

N 

^ = -Y. *^4cj + c]c,) + Y^ e4c, , (2-1) 

(ij) i=^ 

where the first sum runs over the nearest neighbors couples of the lattice, the sec- 
ond sum runs over all N sites and cl,Ci are fermionic creation and annihilation 
operators. For simplicity we consider spinless fermions. The on-site energies ej are 
independent and identically distributed random variables, that we shall take in the 
following uniformly distributed in the interval [—W/2,W/2]. As the fermions have 
no interactions, the study of this model amounts to determine the spectrum of the 
N X N random matrix Hij such that Hij = —t if i ^ j are neighbors on the lattice, 
Hii = Ci on the diagonal. In the following we will take t = 1, which means that e^ 
are measured in units of t. 

As anticipated in the introduction, the lattice we will focus on is a Bethe lattice. 
It is possible, and this is the point of view of almost all mathematical works, to 
define it as an infinite regular tree, i.e. a graph without loop where every vertex has 
the same degree (that we shall denote A; + 1 in the following) . It is however necessary 
for the study of thermodynamical properties of statistical mechanics model (e.g. 
to define a free-energy) to consider finite-size versions of the model, with a given 
number N of vertices, before taking the thermodynamic limit A^ — t- oo. There are 
at least two possible ways to define a finite-size Bethe lattice. The first one is to 
consider a finite tree of depth n, that is a tree in which the root vertex (at generation 
0) has k + 1 offsprings, each of the vertices in the generation from 1 to n — 1 has 
k offsprings, and finally the leaves of the n-th generation have no descendent. The 
second way is to consider a random A; + l-regular graph, i.e. a graph chosen uniformly 
at random among all the graphs on N vertices where each of the vertex has degree 
k + 1. The properties of such random graphs have been extensively studied (see 
ref. |19p for a review). It is known in particular that any finite portion of such 
a graph is a tree with a probability going to one as A^ — >• c«. The advantage of 
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the second construction is the absence of any boundary effect: ah sites, for any 

finite A^, have k + 1 neighbors, whereas in the trees of the first construction the 

sites on the boundary, i.e. on the leaves of the n-th generation, are asymptotically 

as numerous as the bulk (volume) sites and have connectivity one, a pathological 

situation compared to finite dimensions. Random regular graphs can be thought of 

as regular trees wrapped onto themselves. This absence of boundary is particularly 

important for frustrated models, as discussed notably in ref. I20p. 

The spectral properties of H are easily studied in terms of the resolvent matrix, 

which is defined as: 

1 

where I is the identity matrix. For example, the density of eigenvalues Aq,, p{E) = 
J2a ^(^ ~ ^a)/^, can be obtained as: 

p( E) = Urn -—ImTr G(E + if]) . (2-3) 

Moreover the localization properties of the eigenvectors of H can be deduced from 
the behavior of lim rjlGaiE + ir/)P (see for instance refs. l2ip . l22p ). 

r;-s>0 

When the non-zero off-diagonal elements of H have the structure of a finite tree 
it is easy to obtain the diagonal elements of the resolvent from 

j&di 

Here di denotes the set of neighbors of i and Gj^i{z) is the Gjj{z) resolvent of a 
modified matrix in which the edge between i and j has been removed. These new 
variables verify the following recursive equations: 

These equations were first obtained in ref. [Tj) using perturbation theory. Another 
simple way to obtain these equations is to consider a Gaussian model with a kernel 
H — zX (we shall come back on this in Sec. [5]). More formally one can use resolvent 
identities to derive them; this last method has the advantage of being also valid 
directly on the infinite Bethe lattice, the matrix H being turned in an operator 
acting on square-normalizable functions on the sites of the infinite tree. Note that 
it has been proven in ref. |23]) that the random regular graph case can be studied 
using these recursion equations in the thermodynamic limit, thanks to the local 
convergence of random regular graphs to trees. 

Let us denote, with a slight abuse of notation, G{z) the diagonal element of the 
resolvent for an arbitrary site in the Bethe lattice (the infinite tree or the thermody- 
namic limit of random regular graphs), and G{z) the similar quantity in the absence 
of one incident edge. These are random variables, because of the randomness in the 
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choice of the Ha = e^. The above written equations imphes that their distribution 
follows from 



e - ^ - E-=^i' G.(z) ' ^^'^ e-z- Eti Gdz) 



G{z) = ^^iTT^T-, Gi^) = ^. ^,. ^ (2-6) 



where = denotes the equality in distribution of two random variables. In the r.h.s. 
the various Gi are i.i.d. copies of G, and e is independently drawn uniformly on 
[—W/2, W/2]. The spectral properties of the model can be obtained from the solution 
of these distributional equations, by computing 

p{E) = lim -ImE[G(£; + ir?)], L{E) = lim -r?E[|G(^ + ir])\^] , (2-7) 

where the expectation E[-] is with respect the distribution of G. The first quantity 
gives the density of states of energy E, while the second one corresponds to the 
average inverse participation ratio of the eigenstates of energy close to E (for a finite 
graph, the inverse participation ratio of the eigenstate a is L" = ^^ iV'f l^i where ipf 
is the component of the normalized a-th eigenvector on site i): it is non-zero only if 
E falls in an interval of energy comprising localized states (pure-point spectrum) and 
vanishes for extended states (absolutely continuous spectrum). Indeed the density 
of states alone is smooth at the mobility edge separating localized and extended 
states^^' and cannot be used to distinguish the two regimes. An alternative and 
numerically more precise method to compute the mobility edge was proposed in 
ref. [7|) and consists in investigating the stability of a real solution of the equations 
(j'2-6|) for real z = E to the insertion of a small imaginary part; this solution is stable 
whenever E belongs to the localized regime. Physically, for a finite and large tree, 
this amounts to couple boundary sites to a thermal bath and test whether the site 
at the center can dissipate energy at very low temperature when the Fermi energy 
is close to E. If E belongs to the localized regime, energy is not transported across 
the tree and no dissipation is possible far from the boundary, i.e. the imaginary part 
of the self energy vanishes for sites very distant from the boundary. 

A numerical procedure to solve the distributional equations (j2-6p was proposed 
in ref. [T]) , and revivified more recently in the context of finite-connectivity mean- field 
disordered systems under the name of population dynamics,'^ also called the pool 
method.'^ The idea is to approximate the distribution of a random variable, say 
G, by the empirical distribution of a sample of a large number Af of representants 
{Gi, . . . ,Gj\f}. Starting from an arbitrary initial condition a sequence of samples 
is produced. The elements of the new sample G'- are generated by, identically and 
independently for each j, selecting k elements uniformly at random from the current 
sample {Gi, . . . ,Gj^}, drawing an energy e uniformly in [—W/2, W/2], and com- 
puting the value of G"-, according to (|2-6|) . Repeating these steps, one reaches a 
sample approximating the fixed point solution of ()2-6p . and then the observables in 
(j2-7p can be obtained by computing empirical averages over the representants of the 
distribution. The numerical accuracy is controlled by the size M of the samples. 
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§3. Phase diagram, density of states and mobility edges 

We have numericahy solved the equations (|2-5|) with the method described above 
for the connectivity k+1 = 3 and using samples with sizes up to J\f = 10 . The results 
are displayed in Fig. [TJ The innermost solid line is the mobility edge separating the 
regimes of extended and localized states. It crosses the E = axis for a critical 
value of the disorder Wc ~ 17.4, above which no extended states exist. This value 
is in agreement with the one found in refs. l7j).l8|). rrT]) . Moreover for small disorder 
the whole interval [— 2vA:, 2vA;] of the density of states of the pure model {W = 0) 
corresponds to extended states, in other words the mobility edge does not enter the 
band of the pure model up to a strictly positive value of W. This fact has indeed 
been proven rigorously .l^'l^ The dashed lines show the numerically determined 
edge for the density of states, that is the limit of the interval for which p{E) of 
Eq. ()2-7p is found to be non-vanishing. The edge of the spectrum has been proved 
rigorously to be equal to £" = ±{2Vk + W/2) (solid lines in Fig. [1]) by ergodicity 
arguments.'!^ ''^''^ The discrepancy between the numerically computed value and 
the true one is due to the extreme smallness of the density of states close to the 
edge. As we shall discuss in the next section, the decrease of p{E) is extremely fast. 
Therefore, one would need extremely large sample sizes M in order to reproduce the 
density of states close to the edge. We have checked that in agreement with this 
interpretation the location of the observed band edge depends on the sample size J\f 
used. The data presented in the figure for the band edge corresponds to J\f = 2 10^. 

The shape of the density of states and its evolution with the disorder is shown in 
Fig. [21 where we plotted p{E) for W = 0.3 and W = 12. We also show the density of 
states obtained by exact diagonalization for matrices of sizes 8192. The agreement 
between the two methods is very good. However, for strong disorder finite size effects 
become important close to the edge, as expected from the discussion above. In 
Fig. [3] we show the inverse participation ratio L" obtained by exact diagonalization 
(we recall that for a finite system L" = ^^ iV'f l^j where ipf is the value of the 
normalized a-th eigenvector at site i). The evolution of the inverse participation 
ratio with the system size points to a mobility edge £'c — 4 for W = 12, which is in 
agreement with the prediction of the cavity method, see Fig. [TJ 

We have also studied the evolution of these results with the connectivity of the 
Bethe lattice. The phase diagrams are qualitatively similar. The main issue is the 
evolution of the mobility edge with k. In particular we have focused on the critical 
value of the disorder, Wc, above which only localized states exist. This is interesting 
for several reasons. First, it is relevant for recent works studying the effect of electron- 
electron interaction on Anderson localization. In this case one is interested in the 
problem of localization in the Fock space, which has been argued to be related to the 
localization problem on a Bethe lattice with very large connectivity.'^ Second, the 
limit of large connectivity might be the only case where a complete analytical solution 
is within reach. Abou-Chacra et al.!^ developed two different approximation schemes 
to obtain the mobility edges analytically. Both methods lead to an asymptotic form 
of Wc of the form cklnk. The constants c are equal to 4 and 2e ~ 5.43, depending 
on the approximation. We have found that the stability analysis used to determine 
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Fig. 1. Phase diagram for the Bethe lattice with connectivity fc + 1 = 3. The innermost solid 
line indicates the mobility edge between extended and localized states, the outermost solid line 
being the edge of the density of states E — ±(2\/fc + W/2). The dashed line is the numerically 
observed edge, see the text for details. 
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Fig. 2. Density of states, p{E), obtained using the cavity method and exact diagonalization of the 
connectivity matrix of random regular graphs of 8192 sites for W — 0.3 (left panel) and W = 12 
(right panel). The data are averaged over 16 different realizations of the disorder and of the 
graph. 
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Fig. 3. Inverse participation ratio, L, as a function of the energy obtained via exact diagonalization 
for W = 12 and for different sizes of the adjacency matrix (1024,2048,4096). In the locahzed 
regime L does not depend on the size of the matrix, whereas the inverse participation ratio of 
the extended states decrease as the inverse of the system size. The mobihty edge is found in 
i5 ~ 4, in agreement, within the numerical accuracy, with the numerical prediction of the cavity 
method. 



the mobility edge can be interpreted, in the large k limit, as an iterative equation 
on Levy random variables with tail exponent /x = 1/2. This leads to the result 
c = 2\/27r ~ 5.01. Whether this is still an approximation or an exact result is 
unclear at this stage and will be investigated and detailed elsewhere.'^ 

Our numerical results for Wc obtained by solving the distributional equation 
()2-6p are presented in table HI None of the analytical predictions cited above fit 
well the data. Presumably, very strong sub-leading corrections are present, possibly 
scaling like k with potentially In In corrections. In this situation one would need very 
large values of k to find the true asymptotic form. A best fit of our data of the form 
Wc = cklnk-\-bk leads to c = 4.7 and b = 5.9. Imposing one of the three values of c 
quoted above and adjusting the parameter b leads to fits of comparable quality, these 
numerical data are hence insufficient to discriminate between the proposed values of 
c. 



fc + l 


3 


4 


5 


6 


7 


8 


9 


Wc 


17.4 


33.2 


50.1 


67.7 


87.3 


105 


125.2 



Table I. Critical value of the disorder as a function of the connectivity. 
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§4. Lifshits tail and the edge of the spectrum. 

As discussed in the previous section, it was proved rigorously that the edge of 
the density of state is located in E = ±{2Vk + W/2)^'^'^ Numerically we found 
that the density of states becomes extremely small close to the edges. In the following 
we shall develop simple arguments a la Lifshits that allows one to understand the 
result E = ±{2\/k + W/2) and to obtain the form of p{E) close to the edge. 

Let us first recall how the Lifshits' argumenlP^i works on finite dimensional 
lattices. Let us focus on the left edge (the argument for the right edge is the same). 
The key observation is that there is a finite probability to have all the on-site random 
energies inside a sphere of radius R taking values very close to —W/2, say —W/2 + e. 
Using that such a sphere, when disconnected from the rest, is characterized by states 
with energies arbitrarily close to the exact band edge in absence of disorder shifted 
by —W/2 + e, one can construct variational states for the original Hamiltonian that 
have the same energies. This allows one to show that the left and right edges in 
presence of disorder are just the ones without disorder shifted respectively by —W/2 
and W/2, and also to find the decrease of p{E) at the edges. Note that all this can 
be put on a rigorous ground, as done for instance in ref. [5]). 

In this section we shall develop similar arguments to study the behaviour of the 
tail of the spectrum of the Anderson model on Bethe lattices (a related computation 
was presented in ref. [28]) for off-diagonal disorder). Following the strategy developed 
for finite dimensional systems, the first thing to analyze is the eigenvalue problem 
for an infinite tree of connectivity k + 1 without disorder, i.e. where the adjacency 
matrix H is equal to —1 if i and j are nearest neighbours and zero otherwise. We 
look for eigenvectors of H of energy E, i.e., we want to find states {ip) such that 
-H'lV') = E\tjj). We consider spherically symmetric states for which the components 
of I ip) depend only on the distance from the central site of the tree and must satisfy 
the following set of equations: 

Elpn = -ktpn+l - 1pn-l (4-1) 

where Vn is the component of |^) on any given site at distance n from the origin. 
For energies in the range —2\fk < E < 2\/k, the solution of the equations above can 
be written in the form: 

^n = k-l^ f ^ e-^ + i±^ e— ^"l , (4-2) 



with e = {-E + iy/Ak — E'^)/2y/k. The value of c must be chosen in order to satisfy 
the equation for the central site. This yields: 

c = '^-"^ (4.3) 

Finally, the component of the eigenvector 1-0) at distance n from the origin can be 

rewritten as: 

cos(n0) + csinfn^) ,. ,, 
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Thus iV'nP decrease with the distance from the central site as /c~", while the number 
of sites at distance n from the origin grows as fc"-. As a consequence, such radial 
(quasi-)eigenvectors are not normalizable: 

(^1^) = 1 + t±l Y^ p|^^|2 = 1 + ^^ Y^ I cos(n0) + csin(ne)|2 , (4-5) 

n>l n>l 

which diverges for an infinite tree. 

One can however define a modified variational wave- function that is normalizable 
and has variational energy equal to E. The idea is the following: define n^{E) as the 
smallest integer n such that V'n = 0, and consider the truncated wave-function ^^ = 
'i/'n if n < n^{E), ip'^^ = for n > ni,{E). This wave- function is not a solution of the 
equations (|4Tp . yet it is normalizable and verifies /J/ | i/\ = E. The determination 
of the radius n-i,{E) is simplified when E tends to one of the band-edge. Let us focus 
on the left band edge (the argument for the right edge is almost identical) and write 

E = — 2-v/A; + S. At the lowest order in 5 one finds that 9 ~ ^/-^ and c ~ T-rT\/ ^■ 

Solving the equation ^„^ = in this limit one obtains n^ ~ -kJ \/k/5. 

Let us now go back to the disordered case, where the on-site energies Ha are in- 
dependent and identically distributed random variables in the interval [— VF/2, W/2]^ 
and let us use the variational states derived above to determine the position of the 
edge of the spectrum in presence of disorder. There is a finite probability that in- 
side a spherical region of radius n^, all the on-site energies e^ are arbitrarily close 
to —W/2. The variational state constructed above has now a variational energy 
arbitrarily close to —2yk — W/2. The same argument for the right edge leads to 
the 2\/k + W/2. As a consequence, the edges of the Anderson model on the Bethe 
lattice are ^{2yk + W/2), in agreement with the rigorous results based on ergodicity 
theorem.EJ,[I3 

One can make this argument more quantitative and estimate the behaviour of 
the density of states around the left edge of the band, in presence of disorder (again 
the argument for the right edge is essentially identical). Indeed, the probability that 
all the energies in the ball of radius n^ are in the interval [—W/2, —W/2 + 5\ scales 

as {5/W)^~^' *, since the number of sites inside the ball is (— ^) k^* ■ Combining 
this relation with the scaling of n^, with 5 obtained above for the ordered model 
yields the following scaling for the density of states 



p{-2\/k-W/2 + 5) ~exp 



^^^^k-^'^''-''\o^{W/5) 



(4-6) 



This function vanishes indeed extremely fast with 5, as mentioned previously when 
discussing the numerical results of SeciSl Note that, strictly speaking, this expression 
is only a lower bound to the density of states. A form similar to Eq. ()4-6p has already 
been found for models with off-diagonal disorder on Bethe lattices.'^''^ 
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§5. Isolated eigenvalue 

The last point we shall address is whether the smallest eigenvalue Ai is at a finite 
distance from the edge of the density of states at the thermodynamic limit. Note that 
one needs an extensive number of eigenvalues to obtain a p{E) different from zero. A 
single eigenvalue gives a contribution to p{E) of the order of 1/A^, which is negligible 
in the thermodynamic limit. As a consequence, it is in principle possible to have 
the smallest eigenvalue separated from all the others by a gap and the left edge of 
p{E) at a finite distance from Ai. Although this might seem counterintuitive at first 
sight, at least based on the intuition developed for finite dimensional system, it is 
what happens in the case without disorder when H is (minus) the adjacency matrix 
of a random regular graplo- It is easy to check in this case that the constant vector 
is an eigenvector of H with eigenvalue — (fc + 1). It appears thus a gap between 
Ai = — (/c + 1) and the second smallest eigenvalue A2, which concentrates around 
—2Vk, the left edge of the band of the Bethe lattice model, as proved in ref . \T7\i . 

A natural question is therefore what happens to the isolated smallest eigenvalue 
when disorder is introduced. It is also interesting to investigate the properties of 
the corresponding eigenvector, which it is completely delocalized without disorder. 
We shall follow an approach first proposed in ref. I3U|) for off-diagonal disorder, and 
begin by rederiving it with a slightly different presentation, related to the discussion 
of ref. [T8]) . Let us consider a Gaussian probability measure on A^ real variables (pi 
defined by the averages 



(•) = —: / d0i . . . d(j)N • exp 



l^^MHi,j-ESij)(l)j 



2 



(5-1) 



where AA is a normalizing factor. It is clear that the integral is convergent only if 
£^ < Ai, and that by symmetry {(pi) = for all i. If however the limit £^ — >■ A^ is 
taken then an infinitesimal linear term in the action above is enough to produce a 
spontaneous magnetization {(pi) 7^ in the direction of the eigenvector associated 
to Ai, i.e. {(pi) is proportional to its i-th component. This corresponds to a con- 
densation transition on the first eigenvector. Using the replica symmetric cavity 
method^^' (which corresponds here to the Gaussian Belief Propagation'^^') one finds 
that computing the averages with respect to this weight, on a given graph, amounts 
to solve the following recursion equations 

/i,^,(,/),) = -^e-5(^--^)'^? II f d(Pypj,^i{(Py)e-''^^"^"^^' , (5-2) 
where Ni^j ensures the normalization and Ui ^j(( Pi) is the probability measure for 



when the link between i and j is remove d**^! . We parametrize this Gaussian 



*•* For this phenomenon, there is a difference between an infinite random tree and a very large 
random regular graph. Only for the latter the smallest eigenvalue is separated from the others by a 
gap. 

**' The replica symmetric cavity method consists in assuming that the (f>ji are not correlated 
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probability as: 

^^^^m = -^e-^^^'^-'^-^"' . (5-3) 

Plugging this expression in ()5-2p we find: 

Gi^j = — s:^ tt2 n ' y^^i ^ ~ ^ Hij'Gf^iVj'^i , (5-4) 

j'eeAi ^^^"^^ 

which are indeed the same as derived in ref. l3U|) . The equations linking the 'messages' 
G on the various edges of the graph are decoupled from the y's, and corresponds to 
the recursion equations for the resolvents stated in (|2-5p . The linear term coefficients 
y are solutions of a homogeneous set of equations, hence invariant with respect to a 
common multiplication. The smallest eigenvalue Ai can thus be found as the largest 
value of E such that the fixed point y = is stable. 

To investigate the behaviour of Ai for large regular graphs one turns these equa- 
tions in distributional equations, 



1 ^ - \ 



iG,y)= I _^ _, ^ ^y.G.y,] . (5-5) 

As explained in Sec.[2]this kind of distributional equation is easily solved numerically 
by representing the distribution of the random variable (G, y) as a sample of couples 
{{Gi,yi), . . . ,{Gj\f,yf>/)}. These are updated according to a generalization of the 
method explained in Sec. [51 where now the update is made on the couple {Gj,yj) 
and not on Gj alone. The stability of the y = fixed point can then be determined 
by monitoring the norm 'YlijV^ after each iteration, and one identifies the location 
of Ai as the largest value for which this norm decreases upon iterating Eq. ()5-5p . 

The numerical results of this procedure are presented in Fig. [J] for the connec- 
tivity A; -|- 1 = 3. The smallest eigenvalue Ai sligthly decreases upon increasing the 
disorder from the value —3 it has at M^ = 0. We find that even in presence of disorder 
the smallest eigenvalue is separated from the others by a gap. Moreover it remains 
delocalized as shown by its participation ratio that can be obtained from the Vi^j 
when E = X^ . There is a critical value of the disorder, Wi, where Ai becomes equal 
to the edge of the density of states. For larger values of W the smallest eigenvalue 
equals the edge of the density of stateO. In Fig. [J] we show the numerical result 
obtained by the cavity method for the isolated eigenvalue in the case A: -|- 1 = 3, 
which crosses the left edge of the band for a value of the disorder Wi ~ 0.36. 

except through their couphng to i. This is justified in the present case and directly leads to the 
equations (I5-2|) . 

*' For W > W\ our method still finds a value of Ai, which is now larger then the left band edge 
and that approaches the mobility edge when increasing the disorder strength. The interpretation 
of the replica symmetry cavity method in this regime is unclear. It might be that between the 
extensive number of localized states, there is one extended state at least for W not much larger 
than Wi and that the replica symmetry cavity method misses all the other localized states because 
looking at the stability with respect to yi-»j is like studying the stability to a delocalized external 
field, which does not couple to localized states. This point certainly deserves further investigation. 
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Fig. 4. Blow-up of the phase diagram of Fig. [T] the supplementary line starting from —3 is the 
isolated eigenvalue Ai. 



§6. Conclusion 

In this work we have revisited the Anderson localization problem on a Bethe 
lattice. We have shown detailed results on the density of states and the mobility 
edge for connectivity equal to 3. We have also discussed the large connectivity limit, 
whose exact solution remains an open question. Finally, we have shown that for 
not too large disorder the smallest eigenvalue is extended and has a gap from the 
subsequent (localised) eigenvalues. This property may play an important role for the 
disordered Bose Hubbard model on Bethe lattice. Indeed, the usual argument used 
to claim that a Bose glass phase always intervenes between the superfluid and the 
Mott phase may fail. In finite dimensional lattices, when the disorder strength W is 
only slightly larger than the gap of the Mott insulator, it was arguecP^''^ that it 
is favorable to insert particles (or vacancies) but that these are not delocalized and 
hence do not lead to superfluidity since the lowest eigenvalue of the corresponding 
Anderson problem is localized. In the Bethe lattice case, depending on the effective 
Anderson model one finds for these quasi-particles, the smallest eigenvalue may 
correspond to a completely delocalized state with a finite gap from the localized 
ones. Thus, the possibility of a direct quantum phase transition from the superfluid 
to the Mott state is not completely excluded on Bethe lattices. 
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